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JENS MARKLOF 

Abstract. We show that eigenfunctions of the Laplacian on certain non-compact 
domains with finite area may localize at infinity—provided there is no extreme level 
clustering—and thus rule out quantum unique ergodicity for such systems. The con¬ 
struction is elementary and based on ‘bouncing ball’ quasimodes whose discrepancy is 
proved to be significantly smaller than the mean level spacing. 


1. Introduction 

Consider a region T> in M 2 with piecewise smooth boundary and finite area. The 
billiard flow on the unit cotangent bundle of T> is defined as the motion along straight 
lines with specular reflections at its boundary dV. The quantum states and energy 
levels of the flow are determined by the eigenvalue problem for the Dirichlet Laplacian, * 1 


(1.1) ( ( f + ^ = ° 

where A = <9 2 + <9 2 . It is well known that the spectrum is discrete. The asymptotic 
distribution of the eigenvalues 

(1.2) 0 < X 1 < A 2 < ... -> oo 


is governed by Weyl’s law (cf. [SOI 2, E, H Ej and references therein) 


(1.3) lim til :X > <X} = Alea(P) . 

A—>oo A 47T 

The mean spacing between consecutive eigenvalues is therefore asymptotically constant. 
We denote by { <fj}j an orthonormal basis of eigenfunctions, and consider the probability 
measure 


( 1 . 4 ) duj = \<pj(x,y)\ 2 dxdy 

associated with the jth eigenstate. One of the central problems in quantum chaos is 
to classify all weak limits of duj as j —> oo. The quantum ergodicity theorem , due to 
Schnirelman, Zelditch and Colin de Verdiere [23 E2, 13 (adapted for billiard flows on 
domains of the above type in |33j). asserts that, if the underlying dynamics is ergodic, 
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1 Our results can easily be adapted to the case of Neumann boundary conditions provided the spec¬ 
trum of the Laplacian is discrete (which, in contrast to Dirichlet conditions, is not generally the case 

for non-compact regions with finite area). 
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FIGURE 3. Leaky polygonal billiard 


there is a subsequence A^, A J2 ,... of full density 2 such that the corresponding eigen¬ 
functions (fj i (i —> cx)) become uniformly distributed on the unit cotangent bundle of 
V. This implies for instance that for any set A C T> with smooth boundary, 


(1.5) 



Area(Al) 

Area(D) 


The proof of this theorem does not indicate whether in fact all eigenfunctions become 
uniformly distributed (a phenomenon called quantum unique ergodicity since there is 
only one possible quantum limit 123 EH]), or if there may exist sparse subsequences 
that have a singular limit, e.g., measures concentrated on periodic orbits of the billiard 
flow. Such exceptional subsequences have been observed in numerical experiments and 
are referred to as scars or bouncing ball modes. Following earlier results for quantum 
maps HU ESI EDI E3, recent seminal contributions on the question of quantum unique 
ergodicity include the work of Faure, Nonnenmacher and De Bievre mm who prove 
the existence of localized eigenstates for quantum cat maps, and Lindenstrauss’ proof 
|2Bj of quantum unique ergodicity in the case of Hecke eigenstates 3 of the Laplacian on 
compact arithmetic hyperbolic surfaces of congruence type. 

In the present paper we show that for certain non-compact domains V C M 2 with 
finite area the sequence of measures d,Uj is not tight, 4 provided there is no extreme 
clustering of eigenvalues. Hence there exist subsequences of eigenstates g) h that leak to 
infinity, and quantum unique ergodicity is not satisfied for such a system. 


2 A subsequence {Aj 4 }i is °f full density if Hindoo #{i : A J4 < A}/#{j : A j < A} = 1. 

3 Hecke eigenstates are simultaneous eigenfunctions of the Laplacian and all Hecke operators. If the 
spectrum of the Laplacian is simple, as conjectured e.g. for the modular surface, any eigenfunction of 
the Laplacian is a Hecke eigenstate. 

4 A sequence of probability measures dvj is tight if for any e > 0 there is a compact domain K, C T> 
such that limsup^^ J v _ lc dvj < e. 
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Let T> be given by 

(1.6) V = {(x,y) £ l 2 : i > 0, 0 < y < f{x)} 

where / : (0, oo) —> (0, cx)) is right-continuous and decreasing to 0 as x — > oo. More 
specifically, we assume that / is constant on the intervals [a.i,a i+ 1 ), i = 1,2,3,.... 
Examples of such domains are displayed in figs. [T] [31 The condition 

OO 

(1.7) < oo, with 8i := /(a*) and i x := a i+ i - a*, 

2—1 

ensures T> has finite area. To illustrate our main result, let us for example choose 
Si = an d £. — ip where a > p > 0 are abribrary fixed constants. Theorem Q] in 

Section |21 implies that there is a constant C > 0 such that (at least) one of the following 
two statements is true: 

□ There is a subsequence of eigenfunctions ipj t (i = 1,2,...) with eigenvalues 
A j i e 7 r 2 f 2 d +<T ) -f [— Ci~ 2p , Ci~ 2p ] and some c > 0 such that for any compact 
1C CV we have 


( 1 . 8 ) 


lim inf 

2—>00 



> c. 


□ The number of eigenvalues A j in the interval 7r 2 f 2 ^ 1+<T ' ) + [—Ci 2p , Ci 2p ] is un¬ 
bounded as i —> oo. 


The first statement implies that eigenfunctions loose a positive proportion of mass. The 
second alternative implies extreme level clustering; this seems unlikely for a generic 
billiard of the above type, but cannot a priori be ruled out. To get a rough idea on 
whether to expect more level clustering than in the case of compact domains 72, we 
show in Section O that the spectral counting function has the asymptotics (Theorem 12} 


(1.9) 

#0 : Al < A} = - kfiV ~ X + ±VX ± ei ±ljj2r Si Vx) + 0 { V A), 

2=1 r =1 ^ ' 

5iy/\>n 

where 

OO 

(1.10) L{ A) = 2 Y, £ i 

2=1 

is an ‘effective length’ of the boundary &D and J\ is the J-Bessel function. The fluc¬ 
tuations are therefore larger than in the compact case, where the error term is of order 
0(V A); cf. Section [21 for a more detailed discussion. 

The proof of Theorem |l] is elementary and based on the construction of ‘bouncing 
ball’ quasimodes HmHEHIISlEliaElESlIIH] (see also Bogomolny and Schmit’s recent 
work on eigenfunctions in pseudo-integrable billiards 0). The non-compactness of the 
domain allows for quasimodes with discrepancy almost as small as where /./ is 

the quasi-eigenvalue. The best rigorous bound for the discrepancy in the compact case 
is 0(1), cf. H2j. 
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Our construction is completely independent on the choice of / on the interval (0, Gq), 
and one may use this additional freedom to tune / on ( 0 , 0 !) in such a way that the 
billiard flow on T> is ergodic. It seems plausible that this is the case if the billiard flow on 
the restricted compact region Z> 0 = {(x,y) € M 2 : 0 < x < a 1; 0 < y < f(x)} is ergodic 
(as in the examples displayed in figs. |T] and |2J), but to the best of my knowledge there 
are no rigorous results in this direction (see however 23 , 23] GE] f° r proofs of ergodicity 
for different classes of non-compact domains). A further interesting class of examples 
are infinite pseudo-integrable billiards (fig. |3J) that are known to be ergodic 5 for almost 
all initial directions m- 


2. Quasimodes 

A function 0 e H^(D) is called a quasimode for —A with quasi-eigenvalue y and 
discrepancy e, if 


/ 21 n I ||(A + /i)0|| < e||0||, 

W&d = 0 ’ 

where || • || denotes the Lr norm. A sequence of quasimodes {00* with quasi-eigenvalues 
Hi is of order s, if 

( 2 - 2 ) \\(A + Hi)A\\ = 0 (n~ s/ 2 )\\^i\\. 

We summarize a few important properties of quasimodes; more details can be found in 

01221 Ell EH]. 

By expanding 0 in an orthonormal basis of eigenfunctions, 0 = JA(0, it is 

easy to see that m implies 


(2.3) I Wh Vi) ! 2 ( A i - V) 2 A e 2 |M| 2 = e 2 I (Vh Vi) | 2 - 

i 3 

Hence \Xj — y\ < e for at least one j, i.e., there is at least one eigenvalue Xj in the 
interval [y — e, y + e]. Consider the larger interval J = [y — be, y + be\, b > 1. We have 

(2.4) £ KV'ft)! 2 < (be)~ 2 Y. I«W)| 2 (V - U 2 < f>- 2 ||^|| 2 . 

J \j(£J 


For a domain A C T> define 


(2.5) 




|0(x, y)\ 2 dx dy. 


5 Since the modulus of the momentum components in both x- and y-directions are constants of mo¬ 
tion, ergodicity is here understood with respect to a two-dimensional submanifold of the unit cotangent 
bundle. 
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Triangle and Cauchy-Schwarz inequality imply 


( 2 . 6 ) 


A < 


< 


< 




a ? e J 


A 


Xj£j 


A 


Ekt^)I 2 , E 




Aw £ J 


\j(zJ 


X^’^j 

\jij 


'E n« 


JIU 


A 7 G J 


Xj^J 


and hence, together with l ira . 
(2.7) 


'E 

A 7 G J 


2 ^ IMU _ b -l 


QUA — 


IMI 


Now suppose that 

for a sequence of quasimodes ifi with quasi-eigenvalue /i, and discrep- 
(2.8) anq/ e* i/ie intervals Ji = [//j — fee*, /z* + bei\ each contain at most k 
eigenvalues X r 

Then, in each interval Ji there is a such that 


(2.9) 


IQ.IL > Tk 




3. Leaky domains 

Let / : (0, oo) —> (0, oo) be a right-continuous function, monotonically decreasing to 
0 on the half-line [ai,oo) (for some a\ > 0), and f f(x)dx < oo. We are interested in 
the domain D = {(x,y) 6 K 2 : i > 0, 0 < y < f(x)}. In the following we will assume 
that / is chosen so that 

/»oo 

(3.1) / f(x)h(n 2 f(x)~ 2 )dx < oo , 

J a\ 

where h : [0, oo) —> [0, oo) is a fixed increasing function bounded by h(x) < yfx. The 
central result is the following. 6 

Theorem 1. For any given decreasing function r : [0, oo) —> (0, oo), and any infinite 
sequence of real numbers 

(3.2) 0 < pi < y 2 < • • • —> oo 
satisfying 

OO 

(3.3) X T 0^) <oo > 

i =1 

6 The notation A <C B for two positive quantities A, B means there is a constant C > 0 such that 
A < CB. We write AxBifd«B«i. 
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there is a domain V of the above type whose Dirichlet Laplacian has an infinite sequence 
of quasimodes with quasi-eigenvalues 

(3.4) Hi,m,n = n 2 /j,i + m 2 ^i, i,m,n G N, 


and 
(3.5) 
so that 


Ktij) 2 


(!) | (A T Ri,m,n}'fii,m 1 n | O || ||; 

(ii) = 0 for i ± %' OT 71 n', 

(iii) < min{0.001, I m - m / | _1 }||V'i,m,n|| ||^,m',n|| form 7 ^ m', 

(iv) for any compact set K, C V, 


(3.6) 




1 


uniformly for all m,n G N as i —>■ oo. 


Remark 1.1. Note that the set : i,m,n 6 N} is a discrete subset of 

mean density 


with 


(3.7) 

where 

(3.8) 


lim 

A —KX) 


#{(z,ro,w) : < A} _ C_ 

A 47t’ 


c=^-yfi= 


< Area (T>). 


This may either be verihed directly, or concluded from the observation (cf. Sections 0] 
and m that can be identified with the spectrum of the Dirichlet Laplacian on 

an infinite union of rectangles T> t with sides 1 ^ 2 , S, = 7T/q 1//2 , and thus total 

area C = JA Area(Dj). In this interpretation, (14.71) represents Weyl’s law m- 

Remark 1.2. If assumption d2HD holds e.g. for the quasimodes 1 , 1 , eqs. (E2D and 
m imply there is an infinite sequence of eigenfunctions , such that for any compact 

a: c n 


(3.9) 


lim inf W^Wv-k > 


1-6 

Vk 


-l 


That is, the eigenstates Pj % loose a positive proportion of mass. It should be stressed 
that we have not ruled out the probably very remote possibility that assumption m 
with Ci = O(rrfii) can never be satisfied for the domains T> considered in the theorem (an 
explicit construction of T> is given in Section ED. It would be interesting to see whether 
(EED can be established at least for generic choices of such D, i.e., generic choices of Si. 
In Section El we will prove an upper bound for the error term in Weyl’s law, which in 
turn yields a rough estimate on possible level clustering. 
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Remark 1.3. For m, n bounded as i —> oo the theorem establishes quasimodes with very 
small discrepancy, 


(3.10) 


|| (^ “ 1 “ l^i,rn,n)^i,rn,n || O 


Kni 


t (r i,m,n , 

Since h and r can be arbitrarily slowly increasing/decreasing functions (respectively), 
this yields quasimodes of order arbitrarily close to 2; cf. example II. ll below. The number 
of such quasimodes with < A, 

A bb (A) = #{(i,m,n) : m,n = 0(1), < A} 

x : p* < A}, 
is determined by the restriction that 


(3.11) 


(3.12) 


r(A)dA^bb(A) < oo. 


Hence the higher the desired accuracy of quasimodes (achieved by choosing a sufficiently 
slowly decreasing r), the thinner the corresponding sequence of quasimodes becomes. 

Remark 1.4. The theorem also implies that there can be sequences of quasimodes of 
order zero that have almost full density. ‘Order zero’ means that 

(3.13) || (A + Hi, m ,n)i>i,m,n\\ = 0(1) ||^i,m,n||, 

i.e., m£i < Ci for some constant C\ > 0. Since in view of (USD there is a constant 
C 2 > 0 such that £//./,; > 0 2 , we have 

Abb (A) = #{(i,m,n) : = n 2 /.i ?; + m 2 & < A, < Ci} 

X_ 

^ Kni) 2 


(3.14) 


># 


,i, m, n 


n 2 < 


Cf .Ci 

c 2 ,m - 6 


Mi< A 


For suitable choices of h and r this quantity can be arbitrarily close to a function x A, 
cf. (13.1 (SI) . On the other hand, it is bounded from below by a/A. This bound is 
attained in the case when 


(3.15) 


\ ^ v / iu r (/ i X ^ 

Z. «x..v> <0 °’ 


%— 1 


h(tii) 


and coincides with the bound for compact domains, cf. [l3j. Note that the heuristic 
approaches in PEH predict a greater number of bouncing ball modes. 

Example 1.1. Take h(x) = x 13 with 0 < f3 < 1/2. For any given infinite sequence of real 
numbers /i,; with 

(3.16) #{j : /I, < A} x A", 

there is a domain T> with f f(x) 1 ~ 2l3 dx < oo, so that the corresponding quasimodes 'ipj 
have order 2 — 2cr, for any fixed a > 2 (a + (3). That is, 

(3-17) || (A + fJti t m, n )A,m,n II = 0(mfJL^^) HV’t.m.n II, 
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The fact that (f8.1(il) implies (D with t(x ) — x (a' > a) is seen by summation by 
parts. In view of Weyl’s law m and the small discrepancy 0 (/q //}/) for bounded m, 
a failure of assumption (Q would imply an extreme clustering of eigenvalues. As we 
shall see in Section 0 the bounds on the error term in Weyl’a law worsen as a —> 0, 
and hence clustering cannot be ruled out. 

An evaluation of the lower bound for the number of order-zero quasimodes in (18.1411 
yields 

(3.18) Abb (A) > A 0 , 

with 6 = max{l + a — 2a' — 2/3, 1/2}. Note that 9 can be arbitrarily close to 1 for 
suitable parameter choices. 

Example 1 . 2 . A second interesting choice that yields a domain T> with exponentially 
narrow cusps, is h(x ) = y/x/ log 7 (l + x) with 7 > 0. For any given infinite sequence of 
real numbers jtq with 

(3.19) #{j:/b<A}xlog a A, 
there is a domain V with / | log f(x)\~ 1 dx < 00 , so that 

(3.20) ||(AT IXi t m t n)^Pi,m,n\\ 0(?77.1og /7,m,n) 111 /* I > 

for any fixed cr < 2(7 — a). Choose here t(x) = log _ “ x with a' > a, and (18.81) can 
again be checked using summation by parts. 

In this case the number of order-zero quasimodes is bounded from below by 

(3.21) A B b(A) ^ v^A. 


4. Proof of Theorem |U 


We begin by constructing accurate quasimodes on the rectangle [a, a + £] x [0, 5] with 
Dirichlet boundary conditions at y = 0,5. Let y G (^“(M) be a mollified characteristic 
function of the interval [0,1]. That is, 0 < x( x ) A 1, x( x ) = 0 for x [0,1] and x( x ) — 1 
for x G [e, 1 — e] for some fixed, small e > 0. We assume also that x'{ x ) — ^(e^ 1 ) (such 
a choice is always possible). For m, n G N, a e R and i, 5 > 0 put 


^m,n{ x ,y) = X 


(4.1) 
and 

(4.2) = 7 r 2 
Straightforward differentiation yields 


x — a\ . f7rm(x — a)\ . /Trny 


Sill 


Sill 


2 / \ 2 -\ 
m\ f n x 
+ 


(4.3) (A T /lm,n)^m,n(®i 1/) 


1 


2n mx' 


x — a 


cos 


7 T m(x — a) 


i // 

+ X 


x — a 


sm 


7 T m(x — a) 



5 
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and hence 
(4.4) 


|| (A + fim,n)'4^m,n\\ — O ^ 


■m 2 S 


where the implied constant only depends on the choice of x- Because of this and 


(4.5) 

we obtain 

(4.6) 


Nvn|| 2 = -r(l + 0(e)), : 


/ 777- \ 

II 4“ pim,n)'4^m,n\\ O ^ j ||'0m,n|| • 


Furthermore, for n ^ v! we have ('ip m ,n, ipm'y) = 0, and for n = n', m rn !, 
(4.7) 


5 r 

Vpm,ni VW ,n) 7j / X 


2 

5 

2 

£5 


x 


sm 


^ nmx 


r-el /•£ 

+ \ 


'0 


r 1 


X 

XI ~ n 


sm 

2 


nm x 


dx 


. / nmx\ . /Trm'x 


sm 




sm 


dx 


' + 


-fow. 


0 */ 1 —e 


[y(x ) 2 — l][cos(7r(m — m')x) — cos(7r(m + r?7/)x)]<ix 


On the other hand, using integration by parts, we have 
(4.8) f [x(x) 2 — 1] cos( 7 r(m — m')x)dx 

[x( x ) 2 ~ 1] sin(7r(m — m')x) 


7r (m — m!) 


J o 


— / 2x(x)x'(x) sin(7r(m — m')x)dx 


Since y(e ) 2 = 1, sin(0) = 0 the first term vanishes, and since x'i x ) — 0(e 2 ) the integral 
is of 0(1). The analogous argument works for the remaining integrals. Hence 

(4.9) \('>pm,n,'lpm',n)\ < min { 6, - --- 7T l||^ m ,„|| ||^ m ',n||- 

( | m — ntv | J 

We will now give an explicit construction of T>. The function / is chosen constant on 
the intervals [a.;, aj+i), i = 1 , 2, 3,...; set Si = f(a,i) and Pi = a* + i — a*. As quasimodes 
we take 

x — a,:\ . (nm(x - a*) 


(4.10) ipi,m,n( x iy) = X 
with quasi-eigenvalues 

(4.11) = 7T 2 


sm 


m 


sm 


Tiny 


+ 
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By construction, these are completely localized in the rectangle [a;, a*+i] x [0, 5j] and 
hence satisfy requirement (iv) of the theorem. Setting //, = 7 r 2 < 5 F 2 ; every given sequence 
of ^ having property El determines a sequence of Si. Because of El, 


(4 . 12 ) ll(A+ iu’" 1, ’‘ ) |f‘ ,m ’” 11 = °x(™A 2 ) = O x {mStAH) = O x {m^Af). 

To minimize the discrepancy, we would like to choose Ai as large as possible. The choice 
Ai = r(/ij)h(/rj) _1 yields condition (i) and determines /. Since 



f(x)h(n 2 f(x) 2 )dx = ^£ i /(a i )/r(7T 2 /(a i ) 2 ) 

l 


(4.13) 


= T, A,h(^ s i 2 ) 

i 


= < 00, 
i 


the function / is in the required class satisfying El- 
Condition (ii) is evident from El- and (iii) from El- 


5. Asymptotic distribution of eigenvalues 

In view of condition El we would like to control the number of eigenvalues in small 
intervals. The following theorem illustrates that extreme level clustering cannot a priori 
be ruled out. 


Theorem 2. The spectral counting function N( A) = ff{j : A j < A} of the Dirichlet 
Laplacian for the domain V (as in Section w satisfies 

(5.1) JV(A) = X - hhvA + T ItT'-jf irS.Vx) + 0(V A), 

4(/r 47 t 2 tt r V / 

i=1 r= 1 v 7 

Six/~\>7T 


where 

oo 

(5.2) L{ A) = 2 Y, 

i= 1 

5iV\>ir 


and J\ is the J-Bessel function. 


Remark 2.1. The standard bound 


(5.3) 

implies that 

(5.4) 
where 

(5.5) 


J\ (x) -C x 1/2 for x large 

N(\) = Al '° 4 a J P) A + 0 (£(A)vT), 


OO 

i(A) = 2tt J] 


i=l 


1 

V5 


«E 

i=l 

lii<\ 


\[Wi _ 

M/0 
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recall that /q = 7 t 2 /<5 2 and & = 7 t 2 /£ 2 . As the examples following Theorem [T] illustrate, 
a good quasimode discrepancy (£j small) is thus traded with an error bound in (15.411 
approaching o(A). But as we shall see in the following section, cf. eq. (Id.71) . the number 
of eigenvalues in the interval [A, A + a] with a < \/A is 

(5.6) A(A + cr) - A(A) = #{(z, m, n) G N 3 : A < < A + a} + 0 (a/A), 

with quasi-eigenvalues / 4 )TOin as in (E3D- That is, all extreme fluctuations beyond 0( a/A) 
are due to the presence of bouncing ball quasimodes. 


6. Proof of Theorem [2] 

Consider the domains V, — {(x,y) E M 2 : a* < x < a J+ i, 0 < y < f(x)} where 
i = 0,1, 2,... and ao = 0. Let Ay) (A) be the spectral counting function for the Dirichlet 
Laplacian for T> %1 and N^\\) the counting function with Neumann conditions on the 
boundary lines x = a % and x = a l+ \ and Dirichlet conditions on the remaining boundary. 
Set 


( 6 . 1 ) 


n d (\) = Y^n^(\), n k (\) = 


i =0 


i =0 


It is well known (‘Dirichlet-Neumann bracketing’ mm that 

(6.2) A D (A) < A(A) < An (A). 

For i — 0 the general error estimate in Weyl’s law for compact domains yields 

(6.3) Ng\\) = Axe f Do) \ + 0{J\), <>(A)= Am(P 0> A + O(v / A). 

47t 47r 

For the remaining domains we have 

oo 

(6.4) A°(A) := ^a£ } (A) = E N 3 : n 2 y t + m 2 & < A} 

i =1 

and 

OO 

(6.5) A°(A) := A«(A) = A°(A) + #{(n, z) G N 2 : n 2 /q < A}. 

*=i 

Note that 

(6.6) JVj?( A) - JVg(A) <J2\~ = 0 (7A) 
since /.q 1 < cx), cf. (14 . 81) . Therefore 


(6.7) 
Now 

( 6 . 8 ) 


JV(A) = 


Area((Do) 


47T 


A + JVg 1 + 0(V\). 


"S( A)= XI 


i,n=l 


A — n 2 /ij 


& 


+ 0 ( 1 ) 




i,n= 1 
n 2 /ij< A 











12 


JENS MARKLOF 


recall the argument in The main term is 


(6.9) 


oo 

E 

i,n= 1 
n 2 fii<X 


A — 


& 


^EE^ 

»=1 n =1 


n 



OO OO 

= 5^E E 

2—1 n=—oc 

where F((r) = max{l — x 2 , 0}. The Poisson summation formula yields for the sum 

over n 



( 6 . 10 ) 



where F( 0) = 7t/ 2 and for y ^ 0 


( 6 . 11 ) 


So 

( 6 . 12 ) 


F(y) = / Vl — x 2 cos(2nxy)dx 


I -1 


2 y 



The bound (15.31) proves the convergence of the series on the right hand side of (Ena. 
This concludes the proof of Theorem [21 
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